Self-energy-part resummed quark and gluon propagators in a 
spin-polarized quark matter and generalized Boltzmann equations 



A. Niegawa* 

Graduate School of Science, Osaka City University, Sumiyoshi-ku, Osaka 558-8585, JAPAN 

(Received today) 

Abstract 

We construct perturbative frameworks for studying nonequilibrium spin- 
polarized quark matter. We employ the closed-time-path formalism and use 
the gradient approximation in the derivative expansion. After constructing 
self-energy-part resummed quark and gluon propagators, we formulate two 
kind of mutually equivalent perturbative frameworks: The first one is formu- 
lated on the basis of the initial-particle distribution function, and the second 
one is formulated on the basis of "physical" -particle distribution function. In 
the course of construction of the second framework, the generalized Boltzmann 
equations and their relatives directly come out, which describe the evolution of 
the system. The frameworks are relevant to the study of a magnetic character 

of quark matters, e.g., possible quark stars. 
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I. INTRODUCTION 



Recent possible discovery of a quark star [1,2] renews our interest in the study of quark 
matters. It has been pointed out [3] the possibly of existing quark liquid in a ferro-magnetic 
phase. For analyzing the magnetic property of quark matters in a consistent manner [4-6] , 
it is necessary to construct a self-energy-part resummed quark and gluon propagators in a 
spin-polarized quark matter, and thereby frame a perturbation theory. 

The spin-polarized quark matter is, in general, out of equilibrium. For dealing with such 
systems, we employ the closed-time-path formalism [4,5]. In this formalism, propagators, 
vertices, and self-energy parts enjoy (2 x 2)-matrix forms, denoted Let G{x,y) be a 

generic two-point function. Fourier transforming with respect to x — y (Wigner transforma- 
tion), we have G{P, X) with X — {x + y)/2. We assume that G(P, X) depends weakly on X. 
Then, as usual, employing a derivative expansion (DEX), we use the gradient approximation: 

G{P; X) ~ G{P; Y) + {X - Y fdY^G{P] Y) . 

We refer to the first term on the right-had side (RHS) as the leading part (term) and to 
the second term as the gradient part (term). Throughout this paper, we assume that the 
density matrix is color singlet, so that the quark and gluon propagators are diagonal in color 
space and independent of color index. Then, we drop the color index throughout. 

The plan of the paper is as follows. In Sec. II, the leading term in the DEX of the self- 
energy-part resummed (SEPR) quark propagator is constructed. In Sec. Ill, we construct 
the leading term of the SEPR gluon propagator in a Coulomb gauge. In Sees. II and III, 
the argument X is dropped throughout. In Sec. IV, we present the gradient terms of the 
quark and gluon propagators. Then, we frame two mutually equivalent perturbative frame- 
works. One framework is constructed in terms of the "bare" -number-density function (and 
its relative), and the other, which we call physical- A'" scheme, is constructed in terms of the 
"renormalized" -number-density function (and its relative). The latter scheme accompanies 
the generalized Boltzmann equation for the "renormalized" -number-density function and its 
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relatives. The form for the leading part of the SEPR gluon propagator in a covariant gauge 
is given in Appendix D. 



II. QUARK PROPAGATOR 



A. Preliminaries 



Spin-polarization vector 



We define a spin-polarization vector S{P) as follows. For a timehke (P^ — Pq— > 0) 
mode, we choose S^^ = (0, C) (= [(f^ = 1] in the rest frame, where P'^ = (e(po)\/P^, 0). 
Similarly, for a spacelike (P^ < 0) mode, we choose — (0, () in the "po = frame," where 

^ (0, ^/=P2D (= a/=P2^m) [f2 ^ 1^ f (f = 0]. S{P) in any frame, where P'^ = {po,P), 
is obtained through a Lorentz transformation: 



— * 

p- c 


^'^ + e(po)\ 




VP2[\/P2 + 




] 



5''(P) = e{p') { ^ ^ + c'' 



^ , P-C p^ + e(p-ov^^e 



(2.1) 

5-P = 0, = (2.2) 

n" = (1,0). 

When a magnetic field is applied along the ( direction, pq > modes with positive (negative) 
charge go to the state S^{P) (— »S''(P)), while their "antiparticle" counterparts (po < 
modes) go to the state —S^{P) {S^{P)). In what follows, the concrete form (2.1) is not 
used, but only the properties (2.2) will be used. 

The projection operators V±{P) onto the states of definite polarization (±) reads 

^^^p^ ^ l + pe{po)^JiP) 
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Orthogonal basis in Minkowski space and the standard form 



As an orthogonal basis in Minkowski space, we choose 
Pi", SI", 

^ iet'^P^P^NpS^ , (el = -p2Ar2) . 

A generic (4 x 4)-matrix function A{P, N, S, ej_) in a Dirac- matrix space is written in the 
from, 

A = A[ + ^'275 + A'^^ + A'S + A' J + A'^^^ 

+^'775/ + ^^75^ + ^975^ + ^'l075)^± 

+ A'^^fh+ A!^A^$ 
+A\^lfi^^+ A'^^$^^. (2.3) 

We decompose A into four parts, 

^= E ^pAV,= E 'PpA'^'V,, (2.4) 

p, <T=± p, a=± 

and write A'"^ in the form, 

APP = + / + $+Aff$, 

Ap-p = 75 [>ir'' + A^-p t + ^ + ^r'' i^i)^] . 

(2.5) 

It is a straightforward task to obtain 

APP = + pe(po)>l'9 , = ^3 - pe(po)A^'Ai5 , 

= + pe(po)P'^i3 , >if = ^ii + P<V^)A!, , ^2 6) 

AP-P = A'^- pe{po)A', , ^r'' = 4 + P6(po)A;2 , 

APfP = + pe(po)^i4 , APfP = ^lie - P<P^)^io ■ 
We refer Eq. (2.5) to as the standard form (SF) and Ap"^ or AP^ to as an SF-element of A. 
It is to be understood that the (bare and sclf-energy-part rcsummcd) propagators and the 
self-energy part, which appear in the following, are to be written in the SF. 



B. Bare propagator 



First of all, we note that the bare propagator matrix S{P) and the self-energy-part 
resummed propagator matrix G{P) enjoy the "symmetry" property, 



(2.7) 



which results from the hermiticity of the density matrix. Here, fi is the first Pauli matrix, 
t acts on Dirac gamma matrix function, e.g., (A^)^ — A*P'*7j^, and S{P) denotes the 
transpose of the (2 x 2) matrix function S{P), etc. 

The bare propagator S(P) is an inverse of S~^(P) — — 171)%. A general solution to 
^-1^ = SS-^ = 1 is 



S{P)^S^'\P) + Sk{P)M+, 

S^'\P) = E [Sra{P) - MSr - Sa)M^ 
p=± 

Sk{P) = - E Cp-p(P) (Ak(P) - A^(P)) 
p=± 



where the suffix "iiT" stands for the "Keldish component" and 



(2.8) 
(2.9) 

(2.10) 



Sra{P) 



M4 



s 



\ 



R 



\^Sr — Sa —Sa 
(1 ±1^ 



V 



±1 1 

Sr^a) = {t + m)^R^A){P) 



t + m 



p2 _ ^2 _|_ ^PqO+ ' 



fp{p)^e{po)Np{\p^\,p) 



(2.11) 



(2.12) 



Here Sjk^a) is the retarded (advanced) propagator, and iVp(|po|,p) [^p(bo|, — p)] (p = ±) 
is the "bare" number-density function of a quark [an antiquark] with polarization pS{P), 



energy |po| (= + tti^), and momentum p [—p\. Sk in Eq. (2.10) connects opposite 
polarization states. Prom Eqs. (2.7), (2.8), and (2.10), we have 

The derivative expansion is an efficient device for deahng with quasiuniform systems near 
equihbrium or nonequihbrium quasistationary systems. For such systems, Sk is small when 

compared to S^^\ 

C. Dyson equation 

The self-energy-part (E) resummed propagator G obeys the Dyson equation, 

G{P) = S{P)[1 + ±{P)G{P)] = [1 + G{P)±{P)]S{P) . (2.13) 
We write G and E in SFs, 

P, <T=± P, (T=± 

It is worth mentioning that, for the system that enjoys an azimuthal symmetry around 
the (-direction, E, and then also G, are independent of E^, provided that we choose ^ — 
P X C/|P X CI- Then, from Eqs. (2.3) - (2.6), we have 

Er = (p = ±), 

Sr = ^7' (j = 2,3,4). 

Same relations hold for G's. 

Substituting the SFs for S, E, and G in Eq. (2.13), we obtain coupled equations, 

(p,(j-±), (2.14) 

where (stcY" = E^,^=±SP^t^^G'^'' , etc. The relation that involves {..y is to be under- 
stood to hold when sandwiched between projection operators Vp ... Va- We write Eq. (2.14), 
with obvious notation, as 



G = S + SEG = S + GSS , 



(2.15) 



where bold-face letters denote (2 x 2)-matrix in a "polarization space" . 

Prom Eq. (2.7), we obtain the symmetry relations for the SF-elements of (cf. Eqs. 
(2.4) and (2.5)), 



(Gf(P))* = -affiG7(P)fi, 



where 



+ for (P(^,j) = (PP, 1), (PP,2), (pp,3) 
(P - P, 2) , (p - p, 3) , 
(P - P, 4) , 

- for (Pt7, i) = (PP, 4) , (p - p, 1) . 



Similar relations hold for S^^^'s. 



Let us introduce (2 x 2)-matrix function f in the polarization space. 



Then S is written as 



f = diag (/+,/_). 



S = S(°) + SkM+ , 
S(°) = SraI - i{SR - Sa)M+ 



SK = {SR{P)-SA{P))iaC{P), 
C(P) = 



^ C+_(P) ^ 



\^C_+(P) ) 
Among the components of S, is a relation. 



(2.16) 



(2.17) 



(2.18) 



(2.19) 



Then, S is written as 
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S(o) 



± ^ ±(.0) _ ^j^M_ , 



+ — 



= Sii + S21 = — S22 — S12 , 

= fSll - Sllf + fS21 + Si2(l - f) . 

Prom Eq. (2.16) with T,^'^ for G^*^, we obtain the symmetry relations, 
Among the components of G, is a relation, 



Gil + G22 — G12 + G21 . 



Then, G is written as 



G = G(°) + GkM+ , 

Gr 
Gr — Ga —Ga 



V 



(Gflf - iGA)M+ , 



Gr — Gil — G12 — — G22 + G21 , 

Ga — Gil ~ G21 = — G22 + G12 , 

Gk = Giif - fGii + fG2i + Gi2(l - f) 



(2.20) 



(2.21) 
(2.22) 
(2.23) 
(2.24) 



(2.25) 



(2.26) 



(2.27) 
(2.28) 
(2.29) 
(2.30) 



It is worth mentioning that, for equihbrium systems, 'Ek — Gk — 0. Prom Eq. (2.16), 
follows the symmetry relations. 



(2.31) 



Substitution of Eqs. (2.17), (2.20), and (2.26) into Eq. (2.15) yields 

q(0) ^ g(0) _^ s(0)|.(0)q(0) ^ g(0) _^ q(0) j.(0)g(0) ^ 

Gk — Sk + SrYIrGk + Sk'^aGa — SrYIkGa 

= Ok rSk + Gk'^aSa — GrT^kSa ■ 

Prom Eq. (2.32), we obtain 



(2.32) 



(2.33) 



G 



R{A) 



t -m-T. 



R{A) 



(2.34) 



where use has been made of fa. We get from Eq. (2.34), after some 

manipulation. 



-'R 



'^R — '^R ^R ^R — ^R ^R '-^R > 



(2.35) 



(2.36) 



where 



G 



(pre)pp 
R 



m - S^'' 



(2.37) 



As has been remarked above after Eq. (2.14), Eq. (2.37) is to be understood to mean 



P^R 



— m — 



(2.38) 



Such an understanding also applies to Eq. (2.35). Concrete form for Gr^a) will be given in 
the next section. 
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As for Gk, Eq. (2.33), we show in Appendix A that 



Gk = + gP + Gg) , (2.39) 

G^^^ = -Gr1:kGa, (2.40) 
G^K = Gr [757)^C(P)S^ - Sk757/C(P)] G^ 

= GRliiGA, (2.41) 

Gg = Gi^757/C(P) - 75i/C(P)GA . (2.42) 

"Z" of in Eq. (2.41) stands for the "leading part" of the DEX. As mentioned above 
at the end of Subsec. B, for quasiuniform systems near equilibrium or nonequilibrium 
quasistationary systems, G^'* and G^^"* are much smaller than G^'*. The SF for G^'* will 
be given in the next section. The standard forms for H; in Eq. (2.41) and ^^?/[G{P) in 
Eq. (2.42) are also given in the next section. The SFs for G^^^ and G^^ are obtained by 
repeatedly using the formulae in Appendix B. 



D. Self-energy-part resummed propagator G 



It is convenient to introduce 

/ 



G = 



V 



-or 



-ET 



We observe that Eq. (2.35) and G^^"" in Eq. (2.40) are unified to a matrix equation. 



G 



^PP £P-P i (pre) -p-ps-pp 

-m)f3-E -EG E 



(2.43) 



where 



- (pre)pp 

G 



^ ^(pre)pp ^(pre)pp 
Ltr Ui 



'K 



-G 



{pre)pp 



ZPP' 

- m)f3 - E 



(2.44) 
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Forms for G^^l^Jf" and of^^"" 



pp 



The SF for - m)f3 - E reads 



pp 



{t-m)%-t (P) 



~ pp 



mfg + El (P) 



+ 



~pp 

h - ^2 in 



E3 {p)$-t^ {p)f$. 



(2.45) 



One obtains the expressions for the SF-elements G^p^^^ (T — R,A,K and j — 1 — 4) 
through straightforward but tedious manipulation of Eq. (2.44), which includes Eq. (2.38). 
Writing Ej = E^- for short, we have 



G 



G 



G 



(pre)pp 

Rj 

{pre)pp 
Rl 

(pre)pp 

m 



m - 



^ (pre) pp\* 

El 



-ryPP 



G 



1-4), 

(pre)pp _ ^ 
R2 



j)PP 



pre 



G 



(pre)pp 



(Z = 3,4), 



pp 

Kl 



where 



7m{[(m + Ei)2 



(j = l-4), 



pre 



(2.46) 



and 



-£"13 ''-^24 + -E'24''-^13 



2N^Re [FisH24 + F24His] 

— 2P24-RePi^4 , 

2N^F2,ReF[tl , 

—2Re [Fi3 H24 — P24 H13] . 



^13 ^24 — -C'24 ^13 



-2iF2jmF[J^ , 
2zP24/mPS , 



Mi' 



-2F,3ReFl,.l , 
2^iV2Pi3/mP^+3 , 

-E'24 ''-^13 + -E'13 ''-^24 

2iAr2/m[Pi3i/24 + i^24i^l3] , 

-2Pi3PePi423 ) 
2iPi3/mPit34, 

2i/m[-Pi3i/24 + i^24i^l3] , 
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with 



^pre ^pre 

^i?34 = (m + Ei)(l-S;)±7V25.^E:, 
FiSs = (m + Ei)Ea±(l-E2)E*, 

//i3 = (m + Ei)E*, //24 = (l-S2)Ea. 



Form for G"^^^^ and G^^"" , Eqs. (2.35) and (2.40) 



Using the definition (B.l) in Appendix B, one can write the quantity in the square 
brackets in Eq. (2.43) as 



ml fa - E 



■.pp 



t. a (pre) 
E® G 



PP 



(2.47) 



The SF for this is obtained by the repeated use of the formulae in Appendix B: 



Eq. (2.47) = - 
+ 



~ PP ^ ~ —pp ^ ~ —pp ^ ~ —pp ^ ~ —pp' 

~ pp ^ A —pp ^ i —pp ^ —pp ^ A ^PP' 



~pp 



~ -pp 



~ -pp 



~ -pp 



AP-p\^ p2 AP-Pv' ap-py^ JD^ AP-PY" 

3 2 4 3 1 4 2 



~ -pp' 



^4 ^1 ^4 ^2 ^3 ^ ^3 ^2 ' 4 1 



(2.48) 



where 



A?r' 



Al-P 



zp-p a (pre) -p-p 



ZP-P i (pre) -p-p 



, iP-P a (pre) -p-p 



= ^1 ' ■ ■ ■ + P^Es G2 + A^^Eg G3 - P^A^^E, G 



ZP-P i (pre) -p-p 



- p-p i (prc)-p-p i p-p i (pre)-p-p 

E-^ G2 + Sr, 



^2 



„ i p-p a (pre)-p-p - p-p a (pre)-p-p 



ZP-P i (pre) -p-p 

= E^ G3 



. z p-p i (pre) -p-p - p-p a (pre) -p-p 



+ E3 



2~P-P - (pre) -p-p 

-P E^ Go 



i p-p ~ (pre)-p-p i p-p - (pre)-p-p - p-p a (pre)-p-p - p-p - (pre)-p-p 

E^ G4 + E2 G3 — E3 G2 + E4 G^ 



(2.49) 



Wc observe that Eq. (2.48) with Eq. (2.49) is obtained from Eq. (2.45) through the following 
substitutions (T stands for R,A, or K), 
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rpPP I pirppp I rppp -plrppp p2 A/-2 \rppp . rppp rppp , rppp ] 

133 "T 243 "T 313 ~ 423 — -T v [1 144 "T -t 234 ~ 324 "r 414j ) 



^T2 



^T3 



^T4 



rppp I p2rppp I rppp p2rppp 1 , Ar2 rTipp , rppp rppp , rppp 1 

134 "T 244 "T 314 ^ 424 "r ^ L-' 143 "r 233 323 "r 413j ) 



-214 



rppp I pirppp I rppp p2rppp p2 \rppp . rppp rppp , rppp 1 

-'l31"r-' -'241"'"-' 311 -^421] ^ L-' 142 "t" 232 322 "t" 412j ) 



+ 



rpPP I p2rppp I rppp p2rppp , trppp , rppp rppp , rppp 

-'132"'"-' -'242 "T -' 312 ~ -' 422 "T L-' 141 "T -' 231 ~ -'321 "T -'411 



(2.50) 



Here, for T = i? and A, 



and 



DPP 



^p-p/^ipre)-p-p-^-pp 
■'Ri ^Rj ^Rl 



APP _ ^p-Pr^{P'^e)-p-p^-pp 
^iil — ^Ai ^Aj ^Al > 



respectively, and, ior T — 



T^PP 



^Ri 



^(pTc)-p-p^-pp r^{vre)-p-psp-pp 
^Rj ^Kl ^Kj ^Al 



I v^p-p/-((pre)-p-py,-pp 
"T^ATi *^Aj ^Al ■ 

^Wpp 



Then, the expressions for G^-, C^j, and G)^j (j = 1 — 4) are obtained from those of their 



counterparts, in respective order, G^^j'^^^'', G^Jj'^^^^, and G^^^^'^^ with the above substitutions. 



(pre)pp 



The forms for Hf" in d^^" , Eq. (2.4I), and for j5^CP''{P) m Eq. (2.42) 

The form for H^'^ is obtained by using the formulae in Appendix B: 

75 



Hr 



+ (AT^Eg'a - Ar2E^''4J^ - Eg^i)^ + Eg'^^i)^) 
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'p-p 



(2.51) 



The form for -f5liCP^(P) is given by Eq. (2.51) with 

K^j^Q (j = 2-4), 
(i = l-4). 

The form for QP'P 

Having obtained the expression for G^p, we can get the expression for G^'P from Eqs. 
Eq. (2.36), (2.39) - (2.42), by repeatedly using the formulae in Appendix B. 

III. GLUON PROPAGATOR 
A. Preliminary 

We adopt a Coulomb gauge. The result for a covariant gauge is summarized in Appendix 

D. 

As an orthogonal basis in Minkowski space, we choose 

= P^-ponP = (0,p) , e = (O, C - (C •p)p7p') , 
n'^ = (l,0), E^^^e^^'^f-P^CpTi,. 
These vectors are orthogonal with each other and their norms are 

P'^-p', C^ = -l + (C•pT/p^ 

Incidentally, e^-"'^'^ Pp(cr, e^'^^"' PpUa, and e'^'^^'^Cp^o- are not independent but are constructed 
out of the above four vectors, e.g., e^^^P^PpC^ = {E^n" - n^'E'^), etc. 
We define the projection operators, 

VI^^{P) -r--^-^, (3.2) 

r^iP)-^, (3.3) 

n\p) = ■ (3.4) 
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Although, = 1, we have written exphcitly for later convenience. In the above defi- 
nitions, 'T', 'L', and 'G" stand, in respective order, for transverse, longitudinal, and gauge 
fixing. (Following tradition, we call n^ny jv? in Eq. (3.3) the "longitudinal projection oper- 
ator".) Prom Eqs. (3.1) - (3.4), one can show that 

= K'{E±), = SutEI . (3.5) 

Let A be a generic second-rank tensor in Minkowski space, whose (/xz^) -component is 
(A)'"'' ^ Ai^"". A''^ is decomposed as 

U, V=T,L,G 

= E (r^.AuvVvr, (3.6) 

U, V=T,L,G 

+ArE^c , 

/I Ml' _ /iGG'nM'' 
^GG — ^1 I G 1 

A^j^ = Al'^l^vy + , 
A^^Q = ^^J^P'' + Al^^P"" , 

A^'^^ = -Af ^P^n^ . (3.7) 

From Eq. (3.5), follows {Vy ■ Ayy ■ VyT" = A^jy {U, V = T,L,G). We call Eqs. (3.6) and 
(3.7) the standard form (SF) and refer {U,V ^ T,L,G) or A^^ {U,V ^ T,T',L,G) 
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to as an SF-element of A^^^. It is to be understood that the (bare and self-energy-part 
resummed) propagators and the self-energy part, which appear in the following, are to be 
written in the SF. 



B. Bare propagators 



Bare gluon propagator 



First of all, we note that the bare propagator matrix D(-P) and the self-energy-part 
resummed propagator matrix G(P) enjoy the "symmetry" property. 



The first two equations results from the hermiticity of the density matrix. 
D^''{P) is an inverse of 



(D-^P))^" = - 



A 



T3 



with A a gauge parameter. A general solution to (D ^D) = g^^" is written as 

D = D(°) + D;^M+ , 



V 



Dfl 
Dr - Da -Da 



(3.8) 



(3.9) 



(3.10) 
(3.11) 

(3.12) 



(3.13) 



(3.14) 
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where 



/(P) = 9{po)N{\polp) - 9{-po)[l + N{\pol -p)] , 



1 + 

p2 



L 



(3.15) 



Here is the number density of the transverse gluon and Aij(^) is as in Eq. (2.11). From 
Eqs. (3.8), (3.10), and (3.13), we have 

{cry-cr, {cry^-cr (3.16) 

Note that, for the quasiuniform systems near equilibrium, Cj^"^, Cj^\ and C^'^ are small 
when compared to /. 

Bare ghost propagator 
A bare Fadeev-Popov (FP) ghost propagator D is 

f>-j-,rs. (3.17) 



C. Dyson equation 

Gluon sector 

The self-energy-part (H) resummed propagator G obeys 

G(P) = D(P)-D(P)fl(P)G(P). (3.18) 
From Eq. (3.8), we obtain the symmetry relations, for the SF-elements of G'^", 
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G'^^(P) = G[^(-P) {U,V = T,T',L,G), 



where af^ = aj^ with af^ ^+ (U^T,L, G) and 



MV 



+ for {UV,j) = (TL,!), (TG,!) 
- for {UV,j) = (TT',3), (TL,2), 
(TG,2), 



Similar relations hold for (11 



uv 



s. 



Components of n follow the same relation as (2.19) and then n is written as 



n = fi(°) - UkM_ 



n(o) 



/ 



V 



n 



R 



Uk^{i + /)ni2 - Ml . 



(3.19) 



(3.20) 



(3.21) 
(3.22) 
(3.23) 
(3.24) 



From Eq. (3.19) with lij for Gj, we obtain the symmetry relations {U^ V, = T, T', L, G), 



Ky(^'))*=-arn^?(P), 
n^y(^) = n^?(-P) 

-e(po) [N{\po\,p)-N{\po\,-p}] 

X (ulfi-p) - nl^{-p)) . 



(3.25) 
(3.26) 



Components of G follow the same relation as (2.25) and then G is written as 



G = G(°) + GkM+ , 

( ^ . \ 



(3.27) 



G(o) 



Gr — Ga 




Ga 



+ /(G«-G^)M+, 



Gr — Gil — Gi2 , 

Ga — Gil — G21 , 

G^ = (l + /)Gi2-/G2i. 



(3.28) 
(3.29) 
(3.30) 
(3.31) 



For equilibrium systems, = G^^- = 0. Prom Eq. (3.19), follows the symmetry relations 
([/,l^, = T,T',L,G), 



g^/Ap) = gV1{-p) 

+e(po) [N(\po\,p)-N(\po\,-p)] 
X (Glfi-P) - Glfi-P)) . 

Substitution of Eqs. (3.10), (3.20), and (3.27) into Eq. (3.18) yields 
G(o) ^ D(o) _ D(o)nG(o) ^ 

G^ = D,^ - BrUrGk + D^n^GA - BkUaGa 



Eq. (3.36) is formally solved to give 



G 



R{A} 



D"^ + n 



R(A) 



For later convenience, we rewrite Eq. (3.38) as (cf. Eq. (3.9)) 



(3.32) 
(3.33) 

(3.34) 



(3.35) 



(3.36) 



(3.37) 



(3.38) 
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-1 



R 



A 



(3.39) 
(3.40) 



+PQ{n^'P'' + P^n") . (3.41) 

The SF for Gj^^a) will be given in the next section. 

As for Gk, Eq. (3.37), through similar procedure as in the quark case (cf. Appendix 
A), we obtain 



G 



K 



(2) 
K 



G 



(3) 
K 1 



G^^^^GrUkGa, 



p(2) 



G^^ — GrC — CGa , 



-G 



R 



CUa - UrC 



G 



-GrUiGa , 



(3.42) 



where C is as in Eq. (3.14). SF-elements of is given in the next section. The SF elements 
of Gk is obtained by repeatedly using the formulae in Appendix C. 



Ghost sector 



The self-energy-part (ft) resummed propagator G obeys 



G{P) = b{p) 



l + ft(P)G'(P) 



1 + G{P)Ii{P) 



D. 



(3.43) 



Since, D, Eq. (3.17), is a diagonal (2 x 2)-matrix, IT is also diagonal. Then, from Eq. (3.43), 
G is diagonal also. Among the components of IT (G), there is the same relation as (2.19) 
((2.25)). Then we have 

n = ftfa , b = Gt3 (3.44) 

with n and G real, and 



G(P) = 



1 



p2 - n(P) p2 + n(P) 



(3.45) 
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D. Self-energy-part resummed gluon propagator 



Form for (Gr(^a) 



uv 



Through a Slavnov- Taylor identity, {Gi)ug {U — G,T,L) and {G2)tg (and then also 
{G\)gu and {G2)gt via Eq. (3.19)) are related to the self-energy-part resummed FP-ghost 
propagator G and the FP-ghost "pre self-energy part"^ fl^. 

The Slavnov-Taylor identity reads [7]: 



G^iyP'^ — X 



= A 



T3n^ - 



G 



G. 



(3.46) 



Here G{P) is as in Eq. (3.44) with Eq. (3.45). As in Eq. (3.44), 11^ is diagonal (2 x 2)-matrix, 



(n; = n^) . 



(3.47) 



Substitution of the SF for G'"' into Eq. (3.46) yields 



G%?iP) = G'ifiP) = n(P) - G(P) 



tGG 



p2 



p2 ' 



Gl?{P) = GlfiP) = A-^ {E1II,{P)) G{P) , 



Gl^{P) = Gl^{P) = A-^ (C'^n,(P)) G{P) , 
GW,{P) = G\^{P) = A J-^ {nm,{P) -po) G{P) , 



nGG _ r^TG _ r^TG _ r^LG _ n 

'^Kl — '^Kl — '^K2 — '^Kl — U • 



(3.48) 

(3.49) 
(3.50) 

(3.51) 
(3.52) 



is evaluated by replacing the vertex factor gCabcP^ at the "end vertex" with gCabc- Here 
the "end vertex" is the vertex from which the outgoing ghost comes out of the diagram. Then, the 
ghost self-energy part ft is related to fi^ through ft = P^Il^j,. 
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All the above quantities are real. In deriving Eq. (3.48), Eq. (3.45) has been used. Sub- 
stituting Eqs. (3.48) - (3.52) into Eq. (3.35), we obtain = Gf{ = 0. Prom the above 
formulae, we see that Gug (^ — ^> ^) {U = T, L) vanish in the strict Coulomb 

gauge (A = 0). 

We are now in a position to obtain {G^)jjy {U ^ G,V ^ G) from Eq. (3.39). We divide 
{^r)uv ^^^^ pieces, {Guj^y = (G^~°'*^^^ + (^G^''^^^, the latter of which vanishes in 
the strict Coulomb gauge (A = 0). 

Straightforward manipulation of Eq. (3.39) using the formulae in Appendix C yields, for 
the SF-elements of 0^^=°) (= GiJ=°^) (H^^ = H^J), 



^^(A=0)TT 

p^(A=0)LL 
^^(A=0)TL 



VG\ 



{X=0)LT 



(p2 - uT - C'nr) - n^) + P^c^n^uJ^'ul'^ , 
(p2 - n;^^) nf^ - p^c'^'n^^nf ^ , 



(3.53) 



where 



(p2 - n;^^) (p2 - nf^) - p^c'n^n^^n^^] (p^ - n;^^ - c'n^) 



LT 
2 



Here, we note that, from Eq. (3.41), P^ - n'^^^ = P^ - nf^ holds. The SF-elements of the 
gauge-parameter dependent part G^^^ = G^^ reads 

^^(A)TT ^ _p4^~2^2 J|^p2 _ ^'TT ^2^TT\^ | ^p2 _ ^^^^ nf ^ + n^n'^^^nf } 

+C'n2^^ { (p2 _ nf ^) nr + n^nf ^nf } 
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VG\ 



{\)LL 



WLT 



-n^^ {(p2 _ n'TTj (p2 _ n^i^) _ p2^2^4nTLnLr| 

+C'n^^ {(p2 - n;^^) nf^ + p^c^n^n^^n^} 

p2 Lnp (p2 _ n;^^)' + P'c'n'uf^ {n^nr - n^nr^} 



Gf« 



(3.54) 



Here, Gf^ (= G^f) and Gf^ (= G^?) are as in Eqs. (3.49) and (3.51), respectively. G^^ 
(= G%^) is obtained from Gjf , Eq. (3.50), with the help of Eq. (3.32). 



[G^)^y is obtained from th above formulae with the substitutions U^^ 



• Gjj3 , Gjji, and G^2 obtained, in respective order, from G^g , G^| , and G^2 with 
the help of Eq. (3.32) or Eq. (3.34). 



Form for (-^i)^^ i'^ Eq. (3.42) 



Straightforward computation using the formulae in Appendix C yields, for the SF- 
elements of Hi, 
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{m)^ = 2^Im{-c^{Illl+~eIl^s) 

(^or = c'criii'^ . (3.55) 

(^i)f ^, (Hi)^'^, and (^z)^^ are obtained using Eqs. (3.32) and (3.35). Other SF-elements 
than the above ones vanish. 



IV. GRADIENT PARTS OF THE PROPAGATORS AND THE GENERALIZED 

BOLTZMANN EQUATIONS 

Here, we deduce the gradient terms of the quark and gluon propagators, and derive 
generahzed Boltzmann equations and their relatives. Procedure goes parallel to those in 
[8-10], and then we describe briefly. 



A. Quark sector 

A configuration-space counterpart of F{P,X) is denoted by F_{x,y): 

Z(.,.)^/^e-(-)F(P,X) (X^^), 
= {F{P,X)\^^{x,y). 

If F{P,X) is independent of X, F{x,y) = F{x - y) = (F(P))jpT (x - y). Here "IWT" 
("IFT") stands for an inverse Wigner (Fourier) transform. 



1. Preliminary 

Configuration-space counterparts of Eqs. (2.4) and (2.5) are, with obvious notation. 
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P, o-=± 

= + ^ (Af ■ii0) + {i$ ) ■ ) 



(4.1) 



4"-" = 75 



+-(A^-''-(i0)-i+(z,^)-i-^r'') 



Here we have used the short-hand notation F_ ■ G, which is a function whose "(x, y)- 
component" is 



= j d'zF{x,z)G{z,y). 
For later use, we display the Wigner transform of Fjq — F- G: 



Fg{P,X) = F{P,X)G{P,X) - -{F{P,X), G{P,X)} , 



(4.2) 



(4.3) 



which is valid to the gradient approximation. The "Poisson bracket" in Eq. (4.3) is defined 
as 

OF dG OF dG 



{F, G} 



dX^' dP^ dP^ 



(4.4) 



2. Bare propagator and counter-Lagrangian 

We proceed as in [9]. We start from the expression for the free propagator S_{x,y) (cf. 
Eqs. (2.8) - (2.10)), 

S{x,y) = J2^p-S,-r, + Sj,M^, 
p=± 



Sr 



i2.A J 



25 



For the time being, f^{x,y) and Cp_p{x,y) in the above equations are left to be arbitrary. 
Specification of them will be made in Subsubsec. 5. S in Eq. (2.8) is the leading part of the 
derivative expansion (DEX) of S{P,X) S_{x,yj^. Straightforward calculation within 

the gradient approximation yields 



S -S^S-S =1, 

S \x,y) = {i0 ^ - rnj S'^ix - y)f3 - L^{x,y)M_ , 

Llc — LlcI + Llc2 : 



p=± 



+ 



-^c2 = « E ^P75 



p2 ^2 

2 



i E ^P^5 

p=± 



+ 



(4.5) 



(4.6) 



(4.7) 



Here l/P^ = p/p2^ 1/Ar2 = P/Ar2^ and l/e^ = P/ei, with P denoting to take the principal 
part. 

Eq. (4.5) tells us that the free action of the theory [6,9] is 



A = y d^xd'^yip{x)S ^{x,y)ip{y), 



(4.8) 



^ / ^ 
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Since the term with L^(x,y) S. j in is absent in the original action, we should 
introduce a counteraction to compensate it: 

A = / d^xd^y^{x)L,{x,y)Mj{y) , (4.9) 

which yields a (two-point) vertex factor 

iL,{x,y)M_^i{L,,{P,X) + L,2{P,X))M_. (4.10) 

Here "WT" indicates to take Wigner transformation. 



3. Dyson equation 



Let us start with considering a "product" of A and B_ of the type (4.1) (cf. Eq. (4.2)), 



C{^, y)^[A- B] {x, y) 



(4.11) 



Using Eq. (4.3), we obtain, for the Wigner transform of C to the gradient approximation, 

C{P,X)= [VpAP^{P,X)B^'^{P,X)V, 



+^rAA^^^ 



dP^ dXt" 



+ ... 



(4.12) 



where " ... " stands for other pieces of the gradient terms than the second term. Thanks to 
the relation. 



'^^'dP,. dP./^'- 



the second term vanishes, Vp{...}Va = 0. Then, to the gradient approximation, C_{x,y) in 
Eq. (4.11) may be written as 
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C{x,y) ~ 



p,€,,a=± 



Thus, as in Eq. (2.15), we can use the (2 x 2) matrix notation in a "polarization space", 

(Ar=A''-(p,<7 = ±). 

The self-energy-part (^{x, y)^ resummed propagator G obeys 



G = S + S-S-G = S + G-S-S, 

1 



Sr 



J 



S^ = 5^-75l-C-C-75l-5^, 
C 



' c\ 



c_ 



+ 







For G and E, we have (cf. Eqs. (2.26) - (2.30), (2.20) - (2.24)) 



G 



V 



Q-R ~ Q-A ~Q-A 



\ 



[GR-i-i-G^-Gj,]M+, 







y ~^_R + —^A j 

G_K — Q.R • f — f • Q.A + Q-12 ) 
— ~^R • f + f • ^A + ^12 ■ 



Eq. (4.13) may be solved to give 



GR{A){x,y) = l{i0x-m) 6 [x-y) 



-SB(A)(2^,y) 



p [2] , p [3] 



C < [1] 

K — 



(4.13) 



(4.14) 



(4.15) 



(4.16) 
(4.17) 
(4.18) 
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r-[2] 



G 



R 



p[3] 



-SK-C-75iiJ -G^ 

Q.R • H • Q.A ) 

GK-75l-C-C-75ii-G^, 



(4.19) 
(4.20) 



The form for the leading part of the DEX of G{X^ P) G{x, y)^ is the G that is deduced 
in Sec. II. 



4- Gradient piece of the self-energy-part resummed propagator 



Form for Gr^a) 



Prom Eq. (4.16), we obtain, for the component G^, 



X, y) — (i0 — S'^{x — y) — S 



pp 

R 



(pre)-p-p 

(pre)pp 
R 

nPP S^P-P /^(Pi'e)-p-p 



v-ip-p ^(.pre; -p-p v-i-pp 



&R''{x,y) = G^r'''-Tl'R''.G-/-^ 



G^r^''{x, y) = - m) 5\x - y) - Tf^^x, y) 



(4.21) 



(4.22) 



(4.23) 



Solving Eqs. (4.21) and (4.22), we write the solutions as G'^ = G^r^"" ^G^r^"" . Here, G^r^"" 
is the leading part of the DEX of G'^r, whose form has been obtained in Sec. II. The gradient 
part G^R^'^ is obtained as 



q{^)pp _ ;_Q 
R " 



(0)pp 
R 



q{&)pp\-^ (^{^)pp\ 



^R l^il 



1 Q\^)^t^ L 

R j ^ R ( 

X^R ■, j '^R 

^(0)pp /^(0)pp\~^\ ^{0)pp /vP-P|/-'(P 
'^R > y^R ) (~'^R \^R I'-^R 



(pre) 



\^7'\Gr 



s^p-P n 
^R ) ^r 



(pre) 



(pre)' 



(0)pp 
R 

ppr<{^)pp 



^R '^R 



r'{0)pp J v^p-p ^(pre) 
''^R \^R 1 '^R 



l^{0)ppspp-p 
^R ^ 



■'R 



^(pre) P f sp-pp 
^R 1 ^ 



■'R 



G^^"", (4.24) 
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_|_/^(pre)pp ]'j]P-P /^(O)-P-Pl 
,(0)pp 



( r<(^^pp s^p-p\ r<^P^^')~P~p ^ f r<io)pp\s^p-p\r<(p^^)~p~p\ 

X'^R 1 ^R J'^R + \^R I'^R / 



_|_y^(0)pp rj]P-P /^(pre)-p-pl 



where 



{^l^|C} 



„ <9C dA ^ dC 



dX^' dP^ dP^ dXi^ 



Then, Gji{x,y) is written as (cf. Eq. (4.1)) 



p, a=± 



(4.25) 



We write the solution to Eq. (4.23) as G^r^'' = G^P^ + G^P^ ^^^"^ . The form for 

^(pre) (0)pp jj^ ^j^g ^^j.^ ^(pre) (l)pp ^.^^^ ^^ 24) with the 

following replacements, 



G 



(0)pa 
R 



G 



(pre)(0)pCT ^{l)po- ^ ^(pre)(l)pCT ^(pre)=Fp=Fp 



R 



R 



0. 



Ga{P-iX) is obtained from Gr{P,X) with S^'s for Sij's. 



Form for Tlx, which is involved in G^' in Eq. (4- IS) 



Computation of Eq. (4.15) to the gradient approximation yields 



— LlcI + Llc2 ~ X] 2p ■ (^/r/o-)lWT " ifp^ 
P, o-=± 

^ p,a=± 



pcrx _ sr^pa 

A AWT —12 



(4.26) 
(4.27) 



and in Eq. (4.26) come from Ac in Eq. (4.9) (see Eq. (4.10)). The standard forms 
for eJ^, E (Eq. (4.19)), and ^i^ft • C and C • 75]^ in Eq.(4.20) are given in Appendix E. 
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5. Perturbation theories — Generalized Boltzmann equations and their relatives 



The aim of this subsubsection is to construct perturbation theories. We are employing 
the interaction picture in the sense of [11]. Then, the quark-gluon system of our concern is 
characterized by a density matrix p at an initial time — from which fp{P,X^,X) 
and Cp^p{P,Xf,X) (cf. Eqs. (2.8) - (2.10)) are determined. It should be emphasized that 
there is no information at this stage on how fp{P,X) and Cp-p{P, X) evolve in spacetime. 
Then, in the course of construction of a perturbative framework, certain evolution equations 
that describe the spacetime evolution for fp and Cp^p should be settled. As a matter of fact, 
one can choose any forms for the evolution equations, on the basis of which a perturbative 
framework is constructed. Different frameworks are physically equivalent in the sense that 
they lead to the same result for the physical quantities (see below for more details) . In the 
sequel, we construct two kind of perturbative frameworks by employing two different forms 
for the evolution equations. 

As seen from Eq. (4.14), the propagator G is written in terms of Gr, Ga, and Gk- 
Gr{P,X) [Ga{P, X)] is analytic in an upper [a lower] half complex po-plane. Then, in 
calculating some quantity, the parts of G that are proportional to Gr or to Ga yield well- 
defined contributions. Now, we observe that G^' and Eqs. (4.18) and (4.19), contain 
GrGa- Since Ga is essentially the complex conjugate of Gr, GrGa is disastrously large on 
the energy shells^, po = =t'^±(±Pi -^), on which 



i?e(G^/(P,X))-^ 



= 0. (4.28) 

po = ±i^± (±P, X) 



As a matter of fact, in the narrow-width approximation, Im{G'^)~^ —>■ e(po)0''', G^^G'X 
develops pinch singularities at the energy shells in a po-plane.^ Then, G^' and G^' yield 
diverging contribution. In practice, Im {G^^)~^ (oc g"^) is a small quantity, so that the contri- 



^How to find the solution to Eq. (4.28) is given in Appendix F. 
^This is a characteristic feature of nonequilibrium dynamics [12]. 
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butions, although not divergent, are large, which invalidates the perturbative scheme. These 
large contributions come from the vicinities of the energy shells, on which Re (C^j^)'^ ~ 0. 

Appropriate use of the first and second equalities of Eq. (2.36) together with Eq. (2.31) 
shows that G^^'G^''^'' and G'jj'~''G^^^'' do not yield large contributions. This is because, in 
general, the energy shells of Gjf and G^''"'', and of G^^^^^^^ and G^^^^^~^~^, do not coincide. 
For the case of = f_, however, this is not the case. For G^^G^'"'^'', a large contribution 
emerges from the region ReiC^y^ ~ 0, and, for C^^G^"^", large contributions emerge 



(pre)pp^ " 



0. 



from the regions Re (G^'') ^ ~ and from the region Re (G 

From Eqs. (4.18) and (4.19) with Eq. (4.26), we have for the (pp)- and (p— p)-components 
(p = ±) of S - (= G-^' ■ (gW + Gg ) . G-/) , 



WT 



P-d 



'X 



p2 



+ 



N-d 



pe(Po 



Ar2 



lex- 



+ (4.29) 



p-p WT 
K 



m „ , pe(po) 



f P<Po) ,dN, m dN'\ , d 

\ m dPa 2m dPa J dX^ 

jj_dlf__d__ pe{po) ^^ dm d P ■ dx 
\2mdPadX^~^ el 



dPa dX^ 



p2 



Cp-p 



P 



^[(l-^)s^2^ + /psg + (/,-/.)E?r] , 



— 



(75ii-G-E^-S^-G-75l 



pa 



WT 



(4.30) 
(4.31) 
(4.32) 



The first term on the R.H.S. of Eq. (4.29) (Eq. (4.30)) comes from the counter Lagrangian 
L,i (^c2), Eq. (4.6) (Eq. (4.7)), in Eq. (4.26). 

For later reference, we note that the physical number densities, N^^\p,X) and 
N^^\p,X), are obtained through computing current density. 



(fix)) = Tr i^{x)Y^{x)p 



^--Tr{-f^[G^,{x,x) + G,,{x,x)]p} 
Similarly, the physical G±zp(P, X), [P^X)^ is obtained from 



(4.33) 
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--Tr {757^ [G2i(x, x) + G,^{x, x)] p) 



(4.34) 



Bare-N scheme 



As has been emphasized at the beginning of this subsubsection, /p(P, X) in Eq. (4.29) 
and Cp_p{P,X) in Eq. (4.30) {X^ < X^) have not been defined so far. For the purpose of 
determining them, we impose here the condition that the counter Lagrangian is absent, 



P ■ dxff^ = N ■ ax/f ) - ex • ax/f ) = , 

[m{N-dx)+pe{po){e^-dx)]Cf:}p 



pe(po)ex^ 



'2'dK 



dC. 



(B) 



P-P 



dX'^ 



c. 



(B) 
p-P 



(4.35) 



, (4.36) 



where we have written Z^^) (C^.j,) for (Cp_p). Then, Eq. (4.29) (Eq. (4.30)), of 
which the first term on the RHS is absent, is to be solved under the given initial data 
fp{P,Xf,X) {Cp-p{P,X^,X)). Eq. (4.35) is a "free Boltzmann equation" and its relatives 
for the "bare" number densities, N^^\po,p,X) = e(po)f^^\P,X) and N^^Wpo\,p,X) = 
l-^(-Po)/f)(po,-p,X), (cf. Eq. (2.12)). 

The physical number densities, which are obtained from {j^{x)) (Eq. (4.33)), and the 
physical C^p, which is obtained from O's (x)) (Eq. (4.34)), are functionals of and C^^: 

f^^'^\p,x)^eipo)N^^''\po,p) 

+e{-po) [l-7V(P^)(|Po|,-p) 



C^^iiP, X) = GpiP, X; [/f )], [C^%]) . 



(B) 



J-'p and Qp here contain large contributions mentioned above. Solving these equations for 
/^^^ and Cp?p, , one obtains 
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/r = /fH^', X; [^S]) , (4.37) 

= C^J^liP: X- [/(P^)], [CS]) . (4.38) 

In the case of scalar theory [8], the physical number density is shown to obey the generalized 
Boltzmann equation. 

Computation of some physical quantity yields the expression -^([/^^^], [C'p^p]), which 
includes large contribution. Substituting the RHS's of Eqs. (4.37) and of (4.38) for, in 
respective order, f'^^ and Cp^), in F, one obtains the expression F'{[f'^^^]^ [C'p-^^]), which 
does not include large contribution. 

The perturbation theory thus constructed is called the "bare-iV scheme" in [8]. 

Physical-N scheme 

Here we aim at constructing a perturbation theory, on the basis of which no large con- 
tribution appear. Then, in such a scheme, there are no large terms in the relations between 
{fi'^^\ Cl^-i) and {fp,Cp-p). This is achieved if the condition, Eq. (4.29) = Eq. (4.30) = 0, 
could be imposed. This is, however, not possible. Nevertheless, it is possible to construct 
the scheme that is free from the large contributions. 

For determining so far arbitrary fp{P,X) and Cp^p{P,X) [X^ < X^), we impose the 
conditions, 

Tr + Q/(P, Xj) (Eq. (4.29)) = Trlld (Eq. (4.29)) 

= Tr^lfi (Eq. (4.29)) = , (4.39) 

TV75 (/ - m) i)^ (Eq. (4.30)) 

= Tr75 (f^ - nc{P,X)) (Eq. (4.30)) 

= Tr-f5^ (Eq. (4.30)) = . (4.40) 
Here, Qf{P,X) and Qc(P, X) are arbitrary functions with the property, 

{Po = ±<^± (±P, ^) , P, X) 
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= Qc{po = ±uj±{±p, X),p, X) 



1/2 



(4.41) 



As has been discussed in [8], this arbitrariness does not matter (see also subsubsection 6 
below). Computation of Eqs. (4.39) and (4.40) yields, in respective order, 



P ■ dxfp = |Tr + nf{P, X)) [V;" + zEg' 



(4.42) 
(4.43) 
(4.44) 

and 



e± • dxfp = -pe(po)TrJ^^ + ^S 



■'B 



dPry dX'^ 



C, 



p-p 



p-p 

B 



(4.45) 



[mN ■ dx + pe(po)e± • dx] Cp-p 
^~Ti^,(t-nc{P,X)) 



' , .dN. d mdN^ d 
pe{po)e^ ^ ^ 



dPc, dX^ 2 dP^ dX' 



(4.46) 



(4.47) 



These equations are the determining equations for fp and Cp_p, which are to be solved under 
the given initial data, fp{P,X^,X) and Cp-p{P,Xf,X), respectively. 

After imposition of Eqs. (4.42) - (4.47), Rp^p - E^", Eqs. (4.29) and (4.30), turns out 

to 



HPP - = ^ (P^ - nft) Tr (f^" + tT.PP) , 



HP-P - Er - ^75 [P' - 

xTr75(f^-'' + zE^-'') 
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(4.48) 



(4.49) 



On the energy shells po — ±uj±, these quantities vanish, since {P^ — Qf(^c)r){P^ + ^f{c)F) — 
P'^(P'^ — f^/(c7))- Then, the above mentioned large contributions, which turn out to be 
diverging contributions in the narrow-width approximation, do not appear. Thus, + 
turns out to be a well-behaved function. As a matter of fact, in the narrow-width 



approximation. 



+ oc , 



(po ^ ±^±) , 

which is a well-defined distribution. In particular, the relations between the physical 
^y^ph)^ C^^^ and {fp, Cp-p) contain no large term: 

/(P'^) = /, + A/, , gS = Cp^p + AGp_, (4.50) 

with A/p and ISCp^p the perturbative corrections. 

Proceeding as in [9], from Eq. (4.42) on the energy shells, one obtains a generahzed 
Boltzmann equation. In fact, the term with F^'' on the RHS of Eq. (4.42) is proportional 
to the net production rate. To avoid complete repetition, we do not reproduce it here. 



6. Discussion 

Here we like to mention a similarity between the two schemes presented here, the bare- 
scheme and the physical- A" scheme, and those in the ultra-violet (UV) renormalization 
scheme in quantum-field theory. For simplicity of presentation, taking a complex-scalar 
theory, we focus on the mass renormalization and do not mention on the coupling constant- 
and wave function-renormalizations. 



Summary of the UV-renormalization theory 

"Bare" UV renormalization scheme: The free Lagrangian density reads Cq ~ 
—(f)\x){d'^ -\- m\)(f){x) with tub the bare mass. Computation of the physical mass Mph 
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yields Mph = Mph(mB), which includes diverging terms. Solving this equation for mg, we 
have iriB — mB(Mph). Perturbative computation of some physical quantity F yields the ex- 
pression F = F{mB)-i which contains, in general, UV- divergences. Substituting the equation 
ruB — mB{Mpii) for ttib in F{mB), one gets ^^(Mph) = F(mB(Mph)), which is free from 
UV-divergence. 

Physical UV renormalization scheme: One introduces new free Lagrangian C'q — 
— -|-m^)0(x) with m the renormalized mass. Then, the counter Lagrangian should 
be introduced, Cc — Cm — C!m = 4>^{x) \rv? — m%] — is determined so that 

the perturbatively computed physical mass Mph is free from the UV-divergence. Thus, no 
diverging term is involved in the relation Mph = Mph(m). However, there is arbitrariness 
in the definition of the finite part of m, which is determined by imposing some condition. 
This arbitrariness is called a "renormalization scheme dependence" (see, e.g., [13]). It is 
well known that, when one computes some physical quantity F up to, say, nth order of 
perturbation theory, the above arbitrariness affects F at the next to the nth order. The 
renormalization scheme, in which m = M^'^^\ is convenient for many cases. 

Summary of the two schemes presented above 

Bare-N scheme: No counter Lagrangian is introduced. Computation of the physical 
number densities (that are related to /^^^^) and C^p, which are the functionals of /^^^ 
and CpJp, include large contributions. Perturbative computation of some quantity yields 
the expression, which is written in terms of /^^^ and C^p and includes large contributions. 
Rewriting it in terms of the physical quantities, f^^^^ and C^p, one obtains the large- 
contribution free form. 

Physical-N scheme: We introduce a counter Lagrangian — L^^ -\- L^^i which is de- 
termined so that the perturbatively computed physical number densities and C^^'^ do not 
contain large contributions. There is arbitrariness in the definition of the "finite parts" of fp 
and Cp^p. The arbitrariness in the choice of the functions Q,f (Eqs. (4.42) and (4.48)) and 
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Qc (Eqs. (4.46) and (4.49)) is this arbitrariness. It is worth mentioning that, if we could 
choose Q/ and so that fp{P,X) = fj;p^\P,X), Eq. (4.42) on the energy-shell turns out 
to a genuine (generalized) Boltzmann equation. In the opposite case, the function that obey 
the generahzed Boltzmann equation is fp and the physical f^^^^ is written as in Eq. (4.50). 

Similar comment to the above one at the end of the Physical UV renormalization scheme 
may be made here. 

Correspondence 

Above observation discloses the correspondence between the two schemes presented here 
and those in the UV- renormalization scheme: 

Bare scheme: 

/ d^x£o{x) ^ Ao in Eq. (4.8) with L^ = 0, 
TUB ^ /f ) and Cf\ , 
Mph ^ /(P*^) and , 

Physical scheme: 
/ rf^x£o(a;) ^ Aq , 
jd^xLc{x) ^ Ac in Eq. (4.9) , 
rri'^ fp and Cp_p , 
Mph ^ /(P'^) and , 
Absence of 

divergence <-> large contribution , 
Arbitrariness in m <-> Arbitrariness in ft f and Qc , 
Scheme with 

m = M(ph) ^ /, = /(Pi^) and C,_, = C^^'j . 
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B. Gluon sector 



1. Preliminary 



Configuration-space counterparts of Eqs. (3.6) and (3.7) are, with obvious notation, 



U,V=T,L,G j=l 



(4.51) 



where Jtt = 4, Jll = Jgg = Jlg = Jgl = 1 and Jtl = Jlt = Jtg = -^gt = 2, and 

2^''-Ar-(2T)p^+r-4r-r 



4 

E 



7^T^ • iT'^ ■ nV 



~ LL 



LL 



E 

i=i 

2 

E 

2 

E 

2 

E 



TL 



n'^^ ■X':'^ -n^ 



~,TG 



TG 



{id'') A^^ {iff") , 

£l-^f«(za^) + f -^^^(ia^), 

iid^)Af^ . El - {fd^M^ . r , 

nMf^(ia") , 
-(i(?^)Af^n% 



with d'^ — d'' — n'^do and 9o = O/OXq. 



2. Bare propagator and counter-Lagrangian 

We proceed as in [10]. We start from the expression for the bare propagator D(x, y) (cf. 
Eqs. (3.10) - (3.13)), 
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D 



D, 



( 







-Ar + Aa Aa 



-[A«-/-/-A^]M+, 



T3 



IFT 



D.LT — ^TL 



1 



pA 



fa - A 



IFT 



-A 



P4 



IFT 



IFT k \ ^ / ) IFT. 

P^JiFT - - lp4JiFT 



f 3 - A 



Ar ■ Q^^ — ■ A.A 



A r^T'T r<T'T 

iAR ■ Ll3 ~ ^3 



c +c 



c 



Ar ■ C_3^ — Cj^ ■ 



D in Eq. (3.10) is the leading part of the DEX of D{P,X) D{x,y)). Calculation 

within the gradient approximation yields 



V A / IFT 



(4.52) 



'Q2 



fa V^''{id)d^ + V^''{id)d^ - do (d^'rf + n'^S") + V^^iid) ^ + dl 



IFT 

(4.53) 



= 2iM_ 



(P ■ df) V!^" + P- dfni'rf 



^p • ~d~f + ^ao/} (n^p^ + w) + p • dcnp, x)l^c 



-P ■ dCr'iP, X)C^El + P ■ dCP iP, X)E^C 



T'Ti 



(4.54) 



Here P ■ df ^ P^df{P, X) /dX^", etc. 

From Eq. (4.52), we see that the free action of the theory is 
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> = (i^, i^) , (4.55) 

where the color index is suppressed. Eq. (4.55) with Eq. (4.53) tells us that there emerges 
a counteraction, 

which yields a (two-point) vertex factor 

iLr(x,t/)^iLr(P,X)M_. 



3. Dyson equation 



As in Sec. Ill, we use the (4 x 4) matrix notation in Minkowski space. The self-energy- 
part (tl{x, yj^ resummed propagator G{x, y) obeys 



D-D n o. 



D 



D^-D^ -D^ 



+ 



D« /-/ D^ + D^ 



For G and H, we have (cf. Eqs. (3.27) -(3.31) and (3.20) - (3.24).) 



/ 



G 



G 



■R 



y Q-R ^ Q-A —Q-A J 



n 



+ 
/ 



V 

+ 



K 



M. 



+ ) 



Ur 



Ur-1-1-Ua-ILk\m^, 

Q.K = ~Q-R ■ / + / ■ Q-A + Q.12 ) 

n^ = n«-/-/-n^ + ni2. 



(4.56) 



(4.57) 
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Eq. (4.56) may be solved to give 



Q-R(A) — 



R(A) ±^R{A) 

r - n -u r P] p [s] 



G 



G« [n« c-c n^] a 

— G^ • H • G^ , 



G[3] ri f~\ C 
K — ■ i=i ^ Li ■ ) 



(4.58) 

(4.59) 

(4.60) 
(4.61) 



The form for the Icadmg part of the DEX of G{X,P), an Wigner transform of G{x,y), is 
the G that is deduced in Sec. III. 



4- Gradient piece of the self-energy-part resummed propagator 

Form for Gr^a) 

We divide the Wigner transform of Hfk^a) ' G^^^) (cf. Eq. (4.58)) into two pieces (cf. 
Eq. (4.51)), 



ILr{A) • Gk(A) ^RiA)iP, X)Gr(A){P, X) + (llR{A)GR(A)y ^ , 



(4.62) 



(n 



R{A) Gij(A) J 



■hjv -h'v' 

E EE 

U,V,V'=T,L,G j=l j'=l 



R{A)j ^UVt^VV -tjVV 
l<-Rj l<-Lj' ^R{A)j>l<'Rj' 



^Uv'^^^R{A)j^'^R0 '^Lj' '^R{A)j'"^Rj' 



'^'^Lj ^^R{A)j'^Rj '^Lj' R{A)j' ^VV 

■K-Rj' 



dP„ 



^UVfrUV ^UV^VV''^^RiA)j' ^'^Rj' 
-K-Lj ^>-R{A)j '<^Rj l<^Lj' 



dP„ 



(4.63) 
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Here, IlR(^A)iP, X) and Gr(^a){P,X) are, in respective order, Ilji(^A)iP) and Gij(^)(P) in Sec 
III. Using Eq. (4.62) in Eq. (4.58), we obtain the solution for Gr(a) (= G^J^^ + G^JJ^^). 



The form for the leading part G^^^^ is given in Sec. III. The gradient part is 



ji^'^^P -d-^ (P'^a^ + d^P") + |a (P^a^^ + d^'P") I 



p(0) 



- (ns(A)G2(^))^ ^ 



Form /or n^, which is involved in G^^ in Eq. (4-59) 

In the following, we restrict ourselves to the strict Coulomb gauge, A = 0, which is a 
physical gauge. Computation of Eq. (4.57) to the gradient approximation yields 



j UV=T,L 



3 



IWT 



+ (ni2)r • s: 



(4.64) 



In G^''^'^, Eq. (4.59), n^, and then V^" in Eq. (4.64), are sandwiched between G^ and 
G^. Then, in the case of A = 0, P^^P^ and (n^P^ + P^'vy) terms in H^'^, Eq. (4.54), do not 
contribute to G^' (cf. Sec. III). 

The Standard forms for the gradient terms, 11^^'''' and 11^^'^'', and H (Eq. (4.60)) are 
given in Appendix G. 

5. Generalized Boltzmann equation and its relatives 

Structure of the theory is fully discussed in Sec. IVA5, so that we restrict ourselves to 
giving a brief description of the physical- A?" scheme only. 

Same reasoning as in Subsec. IVA apphes here: gJ' and G|J, Eqs. (4.59) and (4.60), 
bring about disaster. This disaster would be overcome if the condition 



n^r - H = 



(4.65) 
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could be imposed. This is, however, not possible. Eq. (4.65) may be imposed for the Vt" , 
n^n", C , and C^E\_ components, which read, in respective order, 

2P . df = -i [(1 + /) (n,,)r - /(n,Or] - i^ll - 2Im {C'ElCr'n^^'^) + ^H^'^^^ , 

(4.66) 

2P ■ df = -i \{i + /) (ni2)t^ - /(n2i)f^^ 



2P ■ dcr = -t 



+ (ni2)r-/(n2i)r 



(1)TT 
K2 



-27m 



2P ■ dcr' 



;i + /) (ni2 



,TT' 



/(n2i)r 



in 



{i)Tr' 



+ iH. 



(l)TT' 



(4.67) 



(4.68) 



(4.69) 



Proceeding as in [10], from Eqs. (4.66) and (4.67) on the energy shells, we obtain the 
generalized Boltzmann equations for the transverse and the longitudinal modes, respectively. 
(As a matter of fact, on the energy shells, the first term on the RHS of Eq. (4.66) (Eq. 
(4.67)) is proportional to the net production rates of the transverse (longitudinal) mode.) 
We do not reproduce them here. It should be remarked that, in the case of L mode, Eq. 
(4.67), the "time-derivative term", dof, in the Boltzmann equation comes from 11^^. More 
precisely, the dof term comes from II^^'"^ (with UV = LL), Eq. (G.l) in Appendix G, which 
is in in Eq. (4.67). Eq. (4.68) [(4.69)] determines spacetime evolution of Cj^ [Cj^'] 
along P. An evolution equation for Cj"^ is obtained from Eqs. (4.69) and (3.16). 

One cannot impose Eq. (4.65) for the remaining C'^n'^, n'^C^, E'^n", and n^E\^ compo- 
nents. This is because, for these components, there are no counterpart of /, C^^, Ci^\ 
and Cj'^. For equilibrium systems, these modes are absent. Then, one can expect that, 
for the quasiuniform systems near equilibrium, these modes do not yield disastrously large 
contributions. 

C. Ghost sector 



The self-energy-part ( Il{x, y) ) resummed propagator G(a;, y) obeys 
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G = D-n-G. (4.70) 

D is an inverse Fourier transform of D in Eq. (3.17). As in Eq. (3.44), D, H, and G are 
diagonal (2 x 2)-matrix functions, D = Dfs =, etc. Solving Eq. (4.70), we see that the 
gradient part of G{P,X) vanishes and 

d{p,x) = (g{p,x) 
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APPENDIX A: RESUMMATION OF THE QUARK PROPAGATOR 

Here we derive Eqs. (2.39) - (2.42). Formally solving Eq. (2.33), we obtain 

Gk — —GrSkGa + G^ , 

G'k = GrS^^Sk (1 + ^aGa) , (A.l) 

where use has been made of Eq. (2.32). Since S^^Sk oc (P^ — w?)5{P'^ — m^) — 0, we 
have G'j^ — 0. This means that the piece Sk of the bare propagator disappears through 
resummation, which is unnatural. 

A correct G'^ is obtained by substituting Eq. (2.18) for in Eq. (A.l) as follows: 

G'k = GrS^^ {Sr - Sa) 75^C(P) (1 + T,aGa) 
= 'GRS-R^SR^,lt - GRSl'^^liSA] C(P) 

X (1 + T.aGa) 
= G^75i/C(P) (1 + S^G^) 

-(l + G«SK)75i/C(P)G^, 
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where use has been made of GrSji^ — {1 + GrEr), which follows from Eq. (2.32). This is 
natural in the sense that 



APPENDIX B: "MULTIPLICATIONS" OF THE TWO STANDARD FORMS FOR 

THE QUARK PART 

We define "multiplications" of the functions of the type (2.5) as the products: 

{A BY" = BP" , [A ® BY" = AP-P B'^" . (B.l) 
Straightforward manipulation yields the SFs (cf. Eq. (2.5)) of {A ® BY" , 

{A ® BYi^' = AYB^^P ± P'^A^YBP^P ± N^AYB^^P 

-P^N^AYBfP , 
{A S)f = AYB^^P ± AYBI^P T N^AYBfP 

+N^AYBt'' , 
{A ® BY^" = AYBt" ± P^AYB^P ± AYBP^P 

-p^AYpp^p , 

{A BYa"^^ = AYBf" ± AYBi"^" T AYBt^ 

+AYB^^'' . (B.2) 

[A BYf" {j = 1 - 4) is given by {A ® B)Y in Eq. (B.2) with BJ^'^ for , and [A B]f 
is given by {A ® BYf" in Eq. (B.2) with BJ"" for Sj""". 
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APPENDIX C: "MULTIPLICATION" OF THE TWO STANDARD FORMS FOR 

THE GLUON PART 



We define a "multiplication" of the functions, A^'^ and B^^'^, of the type, (3.6) with (3.7), 
by C^"" = A^'PBp". Straightforward computation yields the SF for C'"'": 

Cr = A^nr + P^C^n'Al'^Br' + P^~Cn^Al'^Bf - P'C^n' Al"" B^^ , 

-P^~Cn^Al'^Bl^' + n^Al'^Bi^ - P^n'^Al'^Bf + PS^^^^^gt 
-P^Al^B^, 

cr = A^B^' + ~CaYbI^' + AY'bI'^ + n^^f ^^2"^^ - P^^^^Pf^ , 
^ = C'n^^f^Pf ^ + P2^2^^^2''^P2^^ + ^Pf ^ - p2^2^[^Pf ^ , 

^GG ^ _P^^2^2^GT^TG _ p2~^2j^GT^TG _ p2^2 j^GL ^LG ^ J^GG ^GG ^ 

^ = ^r^r + C'^r^r + P^C^n^AY'Bl"^ + ^Pf ^ - p2^^«Pf ^ , 
Ci^^ = ^P^ + C'A^'^Bl'^ + P2^2^2^f Pg^'^ + ^Pf ^ - P^^f^P^^^ , 

C2^^ = c'^r^a"^^' + ^ir^r + ^"^i-^ - p'^if ^pr , 

= ^r^P + ''^i^2^''' + n^yl^Pf ^ + «Pf « , 
Cf ^ = ^Pr + ~CA'i^Bl'^' + n^^f ^P2^^ + «Pf ^ , 

^GT ^ p2^2^2^GT^TT' ^ J^GT ^TT ^ (2 j^GT ^TT ^ j^GL ^LT ^ 
C^G ^ ^2j^LT^TG ^ p2^2^2^LT^TG ^ j^LL ^LG ^ j^LG pGG ^ 
^GL ^ p2^2^2^GT^TL ^ ^2^GT^TL ^ j^GL ^LL ^ J^GG ^GL 
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APPENDIX D: GLUON PROPAGATOR IN A COVARIANT GAUGE 



Here we present a "translation table" to get the expressions for the gluon propagator in 
a covariant-gauge from the Coulomb-gauge counterparts given in Sec. III. 

An orthogonal basis in Minkowski space is given by Eq. (3.1) with the replacements^, 



(D.l) 



h'' = n 



p2 



n 



p2 



I ~ 1 ^ \ 



Then, among the projection operators, Eqs. (3.2) - (3.4), and Vq' are replaced as 



p2 



p2 



Vt is the same as in Eq. (3.2). 
Eq. (3.9) is replaced with 



A 



T3, 



which is already in SF. 

The SF-elements of in Eq. (3.13) are replaced by 

DlliP) = 2mCr{P)e{pMP') , 
Dl%{P)^2mCr'{P)e{po)S{P'). 
Dl\{P) = 27riCf^(P)e(po)<^(P^) , 
Dl\{P) = 2mCl\P)e{po)5{P') , 
DY-{P) = (otherwise) 



4t should be noted that = e^'^P^P^Lha = e^'•'P^P^Cpna. 
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In obtaining these, we have used the fact that [Dug] = [Dgu) =0 [U — T, L), which 
is verified from the "bare counterparts" of Eq. (D.2), below. Dj^ in Eq. (3.15) is replaced 
with 



Eq. (3.16) is replaced with 



{cry 



r'T'T 



Eq. (3.17) is replaced by 

b(p) = 



^ A. ^ 



- -A, 



+ 



J 



Introduction of 11'^, Eq. (3.39), is not necessary, Dg ^ = D and 11'^ = Hr, and, Eqs. 
(3.39) - (3.41) are deleted. 

Description after Eq. (3.43) is replaced with the following one: Solving Eq. (3.43), we 
obtain 



G 



Gr — Ga —Ga 



+ 



f{GR - Ga) + Gk 



M. 



+ ■> 



where 



Gr(P) = G*AiP) = = , 

Gk{P) = -GR{P)flK{P)GA{P) . 

Ur^II*^^ fin + ni2 = -622 - 621 , 

n,, = (i + /)nn-/n2i. 



Eq. (3.46) is replaced with 



G^iyP^ — A 



- P^. 



G. 



(D.2) 
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Eq. (3.47) is deleted. 

Eqs. (3.48) - (3.52) are replaced by 

1 



A 



= A- 



P2 + ipoO+ 
1 



"n(P) - p2] Gr{p) = -A 



1 



P2 + ipoO+ ' 



(P2 + ^PoO+)Sl 

= A .^. , (C'^n^.(i')) G^(i') , 



G^f(P) 
G|?(P) 



A 



-A 



(P2 + ipoO+)C' 
1 

(P2 + ?:7;o0+)n2 

p^' f p2 - n,,) (7,, 



Mk^(P)) ^^(P), 
^kGa 



GJc^{P) — A-^ 



(P2 + iO+)(P2 - iO+) ^ P2 - ipoO+ 

EjJ^Kij. Ga 



0, 



(P2 + iO+) (P2 - iO+) P2 - ipoO+ 

C^'Ur^P^Gk C'^'^k^Ga 



(P2 + iO+)(P2 -iO+) P2-ipo0+ 



^^^(P) = A^ 



"Hi,,, P-Gi 



_(P2 + iO+)(P2 -iO+) P2-ipo0+ 

In Eqs. (3.53) - (3.54), the replacements (D.l) and (n.j)uv^ 

are made, and, 

in the formulae in Appendix C, the replacement (D.l) is made. 



APPENDIX E: STANDARD FORMS FOR THE QUANTITIES IN SEC. IV A 



Standard form for 



From Eq. (4.27), we obtain, after some algebra, 

^ p=± L 

+2^Pe j-^ (P-E^', - N-J:Z - p6(po)elE«'J 



+^P<Po) {-^l^'i^2 - ^l^^m + P6(po)iV2P"E^/,^ 
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where 



2 ap„ 



ON. 



^N^' 



pa — pot _j_ 



p2 



(^^~'')" - ^^«7^ + ^^^7^ 0- = 2, 3, 4) , 



etc. Here, {..., ...} is as in Eq. (4.4). 



[2] 

Standard form for H_ in Gjl in Eq. (4-19) 



Straightforward manipulation of Eq. (4.19) yields 



p, (T=± 

where i?/"^ is as in Eq. (2.51). The SF for each term on the RHS of Eq. (E.l) reads 
fin/f fin y-pp f)M f) r 



+ 



) u y-pp 
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if 



p-p 



p2 



pp 



V,^,lfi {C,_„ E/"} = -P, [n^ E:^f } - TV^J^ E^f } 



— ^ ^ 7, 



Aja^-PP \ S^-PP „,('„ \„t^zllJLs^~PP 



dN, 



a 



p2 



(-P"E:^f + N^Y-,^^ + pe(po)elE:^f ) 



Vp. 
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dC. 



dX 



IdN'^ dN 



2 dPr. 



A4 



dN„, 



■ p2 



IdN"^ dN 



p2 



{C_p„ E^-''} i)^75Pp = [n' {C_,„ ES-3"} - TV^f^ {C,.„ Eg^"} 



IdN'^ dN 



2 JDOS^p-p 

i?4 



" p2 



(P"E£" - iV"E^/ + p6(po)elS?,/) 



p ■ 



Standard forms for 'j^M. ■ C and C • -y^M. in G^' in Eq. (4-20) 



Form for 75ii • C in G^' in Eq. (4.20) is given by Eq. (E.l) with 



E% - 



(i = 2-4), 
(j = l-4). 



C • J5II in Eq. (4.20) is given by Eq. (E.l) with 



^pa 
^Rj 

^ppa 



(j=2-4), 
(j = l-4). 
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APPENDIX F: ENERGY SHELLS OF G''^{P,X) 



To find the energy shells of G^, we need (G^(P, X)) ^, the inverse of Gjf(P, X) (cf. 
Eq. 4.28)). To the gradient approximation, we have 



(F.l) 



Here (gJ^''^) "^ is the (ll)-element of Eq. (2.48) and G^J^^^ is as in Eq. (4.24). If we ignore 
the gradient term in Eq. (F.l), the energy shells are obtained through 



Re{Gf'"'{P,X)y' 



Po = ±Jl\±P,X) 



0, 



where V^p is given by Eq. (2.46) with the substitutions (2.50) being made. Then, the true 
energy shells, po — i<^±(='=P, -'^), are obtained from Eq. (F.l), 



± 



dRe(Gf'"'{P,X)) ' 



po = ±ujf\±p,X) 



(a;±(±p,X)-a;f (±p,X) 



Re 



{Gf''{P, X)y' G«''^(P, X) {Gf'^P, X] 



-1 



P0 = ±4"'(±p,X) 



APPENDIX G: STANDARD FORMS FOR THE QUANTITIES IN SEC. IVB 



Standard forms for I^^'^'^ andl^^^'^ in Eq. (4-64) 



From Eq. (4.57) with Eq. (4.64), we obtain 



j UV=T,L 



(E^.d)f + 2i^Re 

c 



p2 



(G.l) 

ni,^ + n--)(c-a)/ 



C ■ P TfTT , TfTG 
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+ 



^K"-nSr)-nSf}(i^x-a)/ 



1 



p 

+4 fef+L^Kr+n^n)(c-5)/ 



1 

P2 



TTV6 , ^ / n^.!? - ^n^^ ) (E^ . a)/ 



p2 



p2 

+ 2in^n''ReU^^dof 



^l^dj - i ( n^f + ^-^ulA (C ■ 9)/ + i-^ul^{E^ . d)f 



1 



p2 



-El (^,.9)/ 



(G.2) 



[...]'s are obtained using Eq. (3.26). 



Standard dorm for H in Eq. (4-60) 



We write H^" = ^(o)'^'^ + ^(i)'^'^^ ^jth the leading term and H^''^'"" the gradient 

term of the DEX of H'^'^{P,X). Straightforward manipulation of Eq. (4.60) yields 



/ ~ \ TT 



IWT 

~ \T'ri 
Hi 



IWT 



(£T)p^+r 



(^0 



~_ \ TT 
2 



c -c 



IWT 









• r - 




IWT ~ 



~_ \ LT 
2 



IWT 



IWT 

■El, 



~_ \ TT 
3 

~ \TL 
2 



J IWT 
1/ 



n 



IWT 



/ ~ \UV 

with as in Eq. (3.55), and 



tTT 



+zc>rpe [- {cr, ^7.}-~e{cr. Tig}+Ei{cr. n^r-n^r} 
-2^cr(c- 9)ns^ + 2c^cr"(p- H)n5r 



+2(c • P) (cj^'(c- H)nsr - cl'^C- H)nsr 
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(^^""(C • d) + 2ReCr{E^ • d)) U 



TT 
Rl 



2 P2 

-2cj^'{p . d)uij + 2(c • p) [-cr^nij + cnc- mil' 

+Cr{C- d)nll^ + 2i{E^ ■ d)Im (C3^^'n^r) + IcnC ■ d)ReIl 



TT 



p2 

^ {c^s"'^ n-} + ^cJ'^(C- d)ii 



TL 
Al 



where A d B = AdB - A d B, and [...]'s are obtained using Eqs. (3.32) and (3.35). 
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